Observation of a Cooperative Radiation Force in the Presence of Disorder 
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Cooperative scattering of light by an extended object such as an atomic ensemble or a dielec- 
tric sphere is fundamentally different from scattering from many pointlike scatterers such as single 
atoms. Homogeneous distributions tend to scatter cooperatively, whereas fluctuations of the density 
distribution increase the disorder and suppress cooperativity. In an atomic cloud, the amount of 
disorder can be tuned via the optical thickness, and its role can be studied via the radiation force ex- 
erted by the light on the atomic cloud. Monitoring cold 87 Rb atoms released from a magneto-optical 
trap, we present the first experimental signatures of radiation force reduction due to cooperative 
scattering. The results are in agreement with an analytical expression interpolating between the 
disorder and the cooperativity-dominated regimes. 

PACS numbers: 42.50.Fx, 42.50.Vk, 42.25.Bs 



Coherence effects in spontaneous radiation and scat- 
tering of light by small and extended clouds of atoms 
has been studied extensively in the past, starting with 
pioneering work by Dicke [1]. In the context of cooper- 
ative effects in the scattering properties of a large cloud 
of atoms [2M3] and localization of light with cold atoms 
[EL @], the interaction of one photon with a cloud of N 
atoms has seen renewed interest. Since the develop- 
ment of Bose-Einstein condensates (BEC) with ultracold 
atoms, light scattering has been shown to be strongly 
modified in quantum degenerate systems of bosons and 
fermions [7Hll|. Using a simple quantum description of 
the photon scattering by N cold atoms, we studied in a 
recent paper [12| the mechanical impact on cold atoms of 
cooperative scattering of light. We showed that the radi- 
ation force acting on a large cloud of atoms can be dras- 
tically reduced, due to both increased forward scattering 
and a reduced total scattering cross section. Monitoring 
the modification of the atomic motion induced by light 
scattering, this study is complementary to the more com- 
monly investigated features of the scattered light [EL E| 
and may provide a new tool for the experimental inves- 
tigation of cooperativity in the interaction between light 
and cold/ultracold matter, including subradiance and lo- 
calization of light by disorder. 

In this letter we report the first observation of coop- 
erative scattering in a large atomic cloud where disorder 
due to density fluctuations tends to suppress cooperativ- 
ity. The measured reduction of the radiation force is in 
very good agreement with the analytical expression de- 
rived in [12j. We discuss the role of the disorder in the 
radiation force before presenting our experimental setup. 
The experimental results will then be compared to the 
theoretical predictions and we discuss further directions 
of this interesting line of research. 

We consider a cloud of N two level atoms (positions 



Tj, transition wavelength A, excited state lifetime 1/r), 
excited by an incident laser propagating along the di- 
rection (Rabi frequency Qq, detuning Aq, wavevec- 
tor ko). We will restrict our analysis to the low inten- 
sity limit, neglecting non linear effects which could arise 
when more than one atom is simultaneously excited. Us- 
ing the Markov approximation, in which the radiation 
escapes from a large atomic cloud in a time much shorter 
than the characteristic decay time, we derived the fol- 
lowing steady-state expression for the average radiation 
force F c acting on the center of mass of large clouds of 
atoms [121 ] : 
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From the structure 



is the single atom radi- 
function Sjv(k) = 



~k SjLi e ^ k k °^' rj , we derived the average value sn = 
(|5V| 2 )<9,0 and the phase function /n = (I^aH 2 cos 0)o,<f» 
where the average is taken over the total solid angle 
of emission of a photon with wavevector k, at an an- 
gle with ko. For smooth density distributions n(r), 
one can compute the structure function by replacing the 
sum with an integral (s n — » Soo, In — » /oo)- We found 
that F c can be influenced by different effects. On one 
side, the finite extent of the atomic cloud can produce 
strong forward oriented scattering, as in the case of Mie 
scattering or, more precisely, Rayleigh-Debye-Gans scat- 
tering [13|. The balance between the momentum of the 
incident and scattered photons and the atoms indicate 
that for forward emission, the net recoil imprinted onto 
the atoms is vanishing, resulting in a reduced radiation 
force (/oo « Sqo). A different contribution to the re- 
duction of the radiation force can be seen in the pref- 
actor of Eq. (pp), which would appear even in a case of 
isotropic scattering (i.e. « 0). The importance of 
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this prefactor can be understood from the cooperative 
coupling of several atoms into the same vacuum mode 
The number of available modes for large spheri- 
cal clouds (with n(r) oc e _r ^ 2(J b)) can be estimated by 
N m oc (kcjR) 2 , resulting in an atom number per mode 
scaling as N/N m oc N/(k<Jn) 2 . This scaling is conve- 
niently related to the on-resonant optical thickness of the 
atomic cloud b = (3A 2 /2tt) / dz n(0,0,z) = 3N/(ka R ) 2 



Using Soo w l/4(ka R ) 2 and w - 2s^ 
predicts for a smooth density distribution Nsqq = 
and thus the following cooperative radiation force: 
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This expression has been derived assuming a smooth den- 
sity distribution and is consistent with a vanishing radia- 
tion force for very large bo. However, Eq.(j2j) also predicts 
a vanishing radiation force for bo — >• 0, where one expects 
single atom physics to become relevant (and F c — » F\). 
This paradox can be solved by realizing that a continu- 
ous density distribution cannot account for small atom 
numbers because it passes by on the dominant role of 
disorder in the low density limit, not untypical for many 
treatments of coherent scattering of light (3-0, 0, E| . Re- 
placing the sum by an integral in the evaluation of the 
structure function eliminates scattering from microscopic 
inhomogeneities in the cloud of atoms, reminiscent of a 
coarse graining when computing effective indices. As in 
the comparison between cooperative scattering and dis- 
order induced localization [5] , one thus needs to estimate 
the fluctuations of s n and /n in Eq. (pp). One can show 
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FIG. 1: (color online) Comparison between analytical expres- 
sions and numerical evaluation for k(jR = 10 with a configura- 
tion average on 10 realizations as a function of atom number, 
(a) Results for sn (blue circles), /n (red squares), and ana- 
lytical expressions 1/N + Soo (blue line), /oo (red line), (b) 
Forces acting on a cloud of atoms with Ao = —100 T: nu- 
merical evaluation (blue circles) of the average cooperative 
force. The full lines indicate (i) the force in the presence of 
isotropic scattering, i.e. assuming f^ = (green dashed line), 
(ii) the force for continuous density distributions without dis- 
order (red dotted line) and (iii) the total force taking into 
account cooperative scattering and disorder (blue line). 



that to very good approximation [see Fig. DJa)], one has 
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The first term in sn describes scattering by a single atom, 
whereas the second term describes interferences between 
scattering from different atoms. For /n we noticed that 
even though its configuration average is very close to , 
the shot to shot fluctuations of /n are bounded by |/at — 
foo\ < 1/N and have a root mean square scaling as <jf oc 
1/N. It is now straightforward to estimate the average 
cooperative radiation force as a function of on-resonant 
optical thickness bo and size ctr: 
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This expression describes the modification of cooperative 
scattering via its impact on the radiation force taking into 
account both, cooperative effects and the role of disorder. 
This expression allows us to make realistic estimates for 
an experimental implementation to measure the cooper- 
ative modification of the radiation force. In Fig. QJb), 
we highlight the various contributions to the reduction 
of the radiation force for extended clouds of cold atoms 
by choosing kciR = 10 and Ao = —100 Y. 

Let us now describe the experimental setup and proto- 
col to measure the cooperative radiation force (see Fig. 
[2]). Different regimes can be investigated, depending on 
the shape, size, optical thickness of the cloud and the 
laser detuning. In order to illustrate the cooperative scat- 
tering, we measure the radiation force by changing the 
optical thickness of the cloud, keeping all other param- 
eters constant. This is done using the following proce- 
dure. We first load a MOT of 8T Rb from a vapor cell. 
All lasers are tuned close to the D2 line of 8T Rb and 
are derived from distributed feedback laser diodes, conve- 
niently amplified with a tapered amplifier and controlled 
via acousto-optical modulators. Large laser beams and 
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FIG. 2: (color online) Experimental scheme to measure the 
cooperative radiation force. Atoms repumped into the F — 
2 state are pushed by a quasiresonant laser beam, whereas 
atoms remaining in F — 1 only experience free fall. 
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ultraviolet LEDs are used to optimize the loading of the 
MOT, where we trap ~ 10 9 atoms in 2 s. We then apply 
a 50 ms temporal dark MOT period where the intensity 
of the repumping laser is reduced by a factor of 10 and 
the detuning of the cooling laser is increased from —4 T to 
— 8T. This allows to compress the cloud and to produce 
a homogeneous Gaussian shaped distribution of atoms, 
mainly in the F = 1 hyperfme ground state, with a tem- 
perature T « 40 /iK. We then switch off all laser beams 
and magnetic field gradient, leaving the atoms in free fall. 
Hyperfme pumping into the F = 1 state is completed 
during a 0.6 ms molasses period without repumper. Be- 
fore the cloud expands significantly, we apply a repump 
laser for 0.5 ms. Changing the intensity of this repump 
laser allows us to prepare a controlled number of atoms 
in the F = 2 hyperfme ground state without changing 
the volume of the cloud. We thus control the optical 
thickness bo of the almost Gaussian cloud with constant 
size (Jr « 650 fim {kcjR « 5 10 3 ). We then apply an hor- 
izontal, circular polarized "pushing" beam, tuned close 
to the F = 2 — )> F' = 3 transition for 0.8 ms. We adjust 
the intensity depending on the detuning Ao in order to 
scatter around 100 photons per atom in the F = 2 state, 
allowing for a net separation of the cloud displaced by 
the pushing beam (see Fig. [3j). 

Atoms in the F = 1 state are not affected by this 
pushing beam. After an additional time of flight (TOF) 
of 13 ms, we detect all atoms by complete repumping 
and fluorescence imaging at 90° from the pushing beam. 
This experimental protocol allows us to distinguish atoms 
which have been exposed to the radiation force (in F = 2) 
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FIG. 3: (color online) Fluorescence imaging (color-coded scale 
normalized to the maximum signal) of atom clouds in the 
F = 1 (left) and in the F = 2 (right) hyperfme ground state. 
The white squares locate the center of the F — 2 cloud which 
is displaced by an amount depending on its optical thickness. 
For this experiment Ao = —4.2 Y and from top to bottom 
b = 5.5,13.5,17.9,19.1. 



from those unaffected by the pushing beam (in F = 1). 
The spatial displacement Ax of the cloud after a fixed 
TOF is proportional to the average radiation force. 

A typical sequence of observed images at Ao = —4.2 T 
for different repump intensities, and thus different optical 
thickness, is shown in Fig. [3j An average of 100 images 
have been taken to improve the signal to noise ratio, pro- 
ducing a slight broadening of the F = 2 cloud due to in- 
tensity fluctuations of the pushing beam. One can clearly 
distinguish the reduction of the spatial displacement, and 
hence of the radiation force, for an increasing number of 
atoms in F = 2. 

A more quantitative comparison to the prediction of 
Eq. for two different detunings is shown in Fig 3J 
where we plot the radiation force, normalized to unity for 
zero optical thickness, for two different values of the laser 
detuning : Ao = — 1.9 T, and Ao = — 4.2 T. An ab intitio 
estimation of the radiation force from the measured val- 
ues of the intensity, detuning and timing of the pushing 
beam yields a satisfactory agreement within 20%. The 
modified radiation force has a Lorentzian dependence on 
the optical thickness bo [see Eq. Q] with a maximum at 
bo = —12. Fitting the experimental data with the laser 
detuning as the only free parameter, we obtain for the 
two curves A^* = -1.4 ± 0.2 T and A^ = -4.3 ± 0.2 I\ 
We have noticed that the good agreement with the the- 
oretical model is less accurate (by up to 20%) if we use 
a linearly polarized pushing beam, because the Raman 
transitions between different Zeeman levels in the case 
of a linearly polarized beam are not included in our 2 
level model. On the other hand, at a scattering rate of 
& 100 kHz, optical pumping into the \F = 2,rriF = 2) 
state by a circularly polarized beam is only weakly af- 
fected by Larmor precession for residual magnetic fields 
below 20 mG. For laser detunings too close to the atomic 
resonance, we observed an important deformation of the 
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FIG. 4: (color online) Experimental data and fits using the 
cooperative radiation force in the presence of disorder for 
Ao = — 1.9 r (red squares) and Ao = — 4.2T (blue circles), 
with typical error bars from shot to shot fluctuations. The 
shadowed area in the inset corresponds to the nonphysical 
region bo < 0. 
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displaced atomic cloud, indicating the role of multiple 
scattering with different importance in the center and 
the edges of the atomic cloud. The above mentioned 
protocol to measure the cooperative radiation force in 
presence of disorder is however restricted to the quasires- 
onant regime. Indeed, for larger laser detunings, optical 
hyperfme pumping puts more stringent constraints on 
the number of photons which can be scattered on the 
cycling F = 2 — >• F' = 3 transition. We have already 
seen hyperfme optical pumping which produces a tail of 
atoms between the initial cloud, and the atoms in the 
F = 2 state which are pushed by the quasiresonant laser 
beam. Different protocols can be used, but for a quan- 
titative measurement of radiation force, it is very con- 
venient to be able to change the number of interacting 
atoms without changing the volume of the atomic cloud. 
In previous experiments, we have observed a reduction of 
the radiation force for increasing atom numbers by using 
a red detuned dipole trap with a detuning in the range 
of 50 — 200 GHz. Even though strong qualitative evi- 
dence for the reduction of the radiation force has been 
observed, a quantitative comparison with our theoretical 
model has been difficult, partly due to the uncontrolled 
shape of the elongated atomic cloud. 

The good agreement between experiments and theory 
might be somewhat surprising, considering that in our 
model we have assumed atoms at rest and that only one 
photon is present in the atomic cloud at the same time. 
The experiments presented in this letter were performed 
at low intensity (saturation parameter s < 10 -2 ). How- 
ever, as we couple to a complex structure (an extended 
cloud of atoms), it would be interesting to see at what 
light intensities deviations from our linear response model 
("one photon") will appear. We obviously expect such 
deviations when entering the superradiant regime, which 
requires a large fraction of the cloud to be in the excited 
state at the same time. Another important assumption 
is that we neglect any Doppler broadening which could 
lead to a dephasing between all induced dipoles [lij . Ex- 
periments on superradiant scattering performed in BECs 
[lOj have for instance noticed that due to such dephasing 
processes no signs of superradiance were observed above 
the transition temperature for BEC. In contrast to our 
experiments, such superradiance with escape rates scal- 
ing as N 2 requires a strong excitation with about half 
of the atoms excited at the same time [1]. Also, in our 
geometry, the cooperative scattering is strongly forward 
directed, reducing the impact of Doppler shift (ko — k) • v 
compared to a geometry where cooperative scattering oc- 
curs at 90° from the incident beam. 

In conclusion, we have derived an analytical expression 
for the radiation force taking into account cooperative 
effects as well as disorder. This allows for a continuous 
description for the radiation force between microscopic, 
single atom Rayleigh scattering and macroscopic cooper- 
ative scattering. We have experimentally observed this 



effect which manifests itself by a reduction of the radi- 
ation force with increasing atom number, well described 
by a simple formula taking account of the finite optical 
thickness of the cloud. The experiment thus represents 
the first measurement of the impact of cooperativity and 
disorder on radiation force. Understanding how to con- 
nect disorder and cooperativity can be important for the 
study of localization of light in cold atoms [5|, |6| . In this 
context the role of increased spatial densities will need to 
be considered, where the phase shift of the driving and 
scattered fields by the atomic cloud will become increas- 
ingly important. We also expect the effects studied in 
this letter to have sizable consequences on the heating 
rate in optical dipole traps. One ingredient in the model 



is the driven timed Dicke state [12| , for which important 
entanglement properties have been predicted [H, Elj . As 
a classical description of many driven dipoles also yields 
strong forward scattering [ljjj], it would be interesting to 
design an entanglement witness [Hj to probe the degree 
of relevant entanglement when N atoms are excited by a 
laser. 
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